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DENSITY PROFILES OF DARK HALOS FROM THEIR MASS ACCRETION 
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M. A. Alvarez, 1 ' 2 K. Aim, 1 and P. R. Shapiro 1 



o 
o 

(N 



> 

m 
m 

(N 
O 

o 

or 

6 
U 

CO 

c3 



13 



RESUMEN 

El resumen sera traducido al espaiiol por los editores. We use the universal mass accretion history 
recently reported for simulations of halo formation in the cold dark matter model (CDM) to analyze the 
formation and growth of a single halo. We derive the time-dependent density profile three different ways, based 
upon three approximations of successively greater realism: equilibrium, radial orbits, and a fluid approximation. 
For the equilibrium model, the density profile is well-fit by either an NFW or Moore profile over a limited range 
of radii and scale factors. For the radial orbit model, we find profiles which are generally steeper than the NFW 
profile, with an inner logarithmic slope approaching -2, consistent with a purely radial collisionless system. In 
the fluid approximation, we find good agreement with the NFW and Moore profiles for radii resolved by N- 
body simulations (r/r2oo > 0.01), and an evolution of concentration parameter nearly identical to that found 
in N-body simulations. The evolving structure of cosmological halos is therefore best understood as the effect 
of a time-varying rate of mass infall on a smoothly distributed, isotropic, collisionless fluid. 

ABSTRACT 

We use the universal mass accretion history recently reported for simulations of halo formation in the cold dark 
matter model (CDM) to analyze the formation and growth of a single halo. We derive the time-dependent 
density profile three different ways, based upon three approximations of successively greater realism: equilib- 
rium, radial orbits, and a fluid approximation. For the equilibrium model, the density profile is well-fit by 
either an NFW or Moore profile over a limited range of radii and scale factors. For the radial orbit model, we 
find profiles which are generally steeper than the NFW profile, with an inner logarithmic slope approaching 
-2, consistent with a purely radial collisionless system. In the fluid approximation, we find good agreement 
with the NFW and Moore profiles for radii resolved by N-body simulations (r/r2oo _ 0.01), and an evolution 
of concentration parameter nearly identical to that found in N-body simulations. The evolving structure of 
cosmological halos is therefore best understood as the effect of a time- varying rate of mass infall on a smoothly 
distributed, isotropic, collisionless fluid. 
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1. INTRODUCTION 

Current understanding of dark matter halos relies 
upon N-body simulations of collisionless cold dark 
matter (CDM) with Gaussian-random-noise initial 
conditions. Two "universal" profiles bracket the re- 
sults ; (Navarro, Frenk, & White 1997; Moore et al. 
1998). The NFW(Moore) profile has an inner den- 
sity profile p oc r -1 (r -1-S ). Wechsler et al. (2002) 
(WBPKD hereafter) found that the mass and con- 
centration of individual N-body CDM halos grow 
over time according to simple universal formulae. In 
what follows, we attempt to explain this result. 



2. THREE DIFFERENT MODELS 

We attempt to understand the form and evo- 
lution of dark matter halos with three spherically- 
symmetric models; equilibrium, radial orbits, and an 
isotropic fluid. Each model assumes that the mass 
M v ir within an overdensity <5 V i r follows the relation 
given by WBPKD 



M vll (a) = Moo exp [-Sa c /a] 



(1) 
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where a c is the scale factor at collapse and S is 
the logarithmic mass accretion rate d(lnAf v ; r ) /d(lna) 
when a — a c . Such a relation is claimed to be a good 
fit to the evolution of halos of different masses and 
formation epochs. As in WBPKD, we have chosen 
5 = 2 and £ V ir = 200, so that the halo has a mass 
M200 and radius r2oo • We have found an initial per- 
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Fig. 1. (Left) Density profile from equilibrium model 
along with best-fitting NFW profile for this profile at 
present. Inset in upper-right shows same over much 
larger range. (Right) Evolution of NFW concentration 
parameter in the equilibrium model. Different line types 
indicate different ranges a?i n < x < 1, within which halo 
was fit to an NFW profile, where x = r/r v i T ,r v i r = r2oo- 



turbation profile consistent with equation (1), 
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where <5; depends on the initial scale factor a; and 
<5 v i r , and M is the unperturbed mass. The free pa- 
rameter b is unity in the absence of pressure or shell 
crossing inside of r v ; r . 

2.1. Equilibrium 

In the simplest model, we have made the assump- 
tion of complete equilibrium inside the halo, so that 
the velocity is zero for r < r v j r = r 2 oo- The mass of 
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□ Radial Orbits : 
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M/M -exp[( 2a c )/a] 




/ □ Fluid Approximation : 
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Fig. 2. Evolution of mass for the radial orbits (top) and 
fluid approximation (bottom) simulations. Shown above 
each are the fractional deviations A = (M oxact —M)/M. 
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Fig. 3. Density profile at the end of the radial orbit 
simulation. 



the halo is given by 



47T 

T 



^vir(a) = -^virPf^ir, 



(3) 



where p is the mean density at that epoch. Mass con- 
tinuity implies the density p v ; r just inside the virial 
radius is 



dM v 



4irp v 



dr v 



(4) 



da da 
Differentiating equation (3) and combining with 
equations (1) and (4) yields 
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where po is the mean background density at a 
The virial radius is given by 



^"vir 



aexp 
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(5) 
= 1. 

(6) 



Equations (5) and (6) are parametric in a, implying a 
radial density profile p(r) = /3 v ir(^vir) which is frozen 
in place as matter crosses r v j r . Taking the limit in 
which a — > oo, the outer density profile approaches 
p oc r _ at late times, consistent with finite mass, 
while the inner slope becomes asymptotically flat. 
The NFW profile is given by 



p{x) 



3x(l + exf 



where 
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(8) 



5(C) In(l + c)-c/(l + c)' 

r/r v j r . Combining equations (5) and (7) 
1, yields an equation for the evolution of 
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Fig. 4. (Top) Density profile at the end of the isotropic 
fluid calculation. (Bottom) Circular velocity profile. 

concentration with scale factor (see Fig. 1), 



2.2. Radial Orbits 

We use a finite-difference spherical mass shell 
code to follow the evolution of a small amplitude per- 
turbation which is chosen so that the resulting virial 
mass will evolve according to equation (1) from an 
initial perturbation given by equation (2). The shell 
code has an inner reflecting core and the results pre- 
sented here used 20,000 shells (see Figs. 2 and 3). 

2.3. Fluid Approximation 

The collisionless Boltzmann equation in spher- 
ical symmetry yields fluid conservation equations 
(7 = 5/3) when random motions are isotropic. Halos 
in N-body simulations have radially biased random 
motion, but the bias is small, especially in the cen- 
ter. This model is therefore a better approximation 
to halo formation in N-body simulations than one 
with purely radial motion. We use a 1-D, spheri- 
cal, Lagrangian hydrodynamics code as in Thoul & 
Weinberg (1995), using 1,000 zones logarithmically 
spaced in mass (see Figs. 2, 4, and 5). The ini- 
tial conditions were chosen in the same way as those 
for the radial orbit model (Eq. 2), with zero initial 
temperature. 

3. RESULTS 

• Equilibrium model docs not reproduce linear 
evolution of concentration parameter with scale 
factor reported by WBPKD, but can be fit by 
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Fig. 5. Evolution of concentration parameter with scale 
factor in the isotropic fluid calculation. 

an NFW profile over a limited range of radii and 
scale factors (Fig. 1). 

• Mass evolution for a perturbation given by equa- 
tion (2) is close to that of equation (1) (Fig. 2). 

• The radial orbit model fails to reproduce the 
inner slope of the NFW profile, approaching 
p oc r~ 2 instead, consistent with the argument 
of Richstonc & Tremaine (1984) (Fig. 3). 

• Fluid approximation halo is well-fitted by NFW 
and Moore profiles for radii resolved by N-body 
simulations (r/r 2 oo > 0.01) (Fig. 4). 

• Evolution of NFW concentration parameter in 
the fluid approximation is a close match to that 
of WBPKD, with c NFW = 4.25a/a c a good fit 
(Fig. 5). WBPKD reported c NFW = 4.1o/a c . 
Complicated merging process is not necessary 
in order to understand density profile evolution. 
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